Quantum information and special relativity 
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Relativistic effects affect nearly all notions of quantum information theory. The vacuum behaves as 
a noisy channel, even if the detectors are perfect. The standard definition of a reduced density matrix 
fails for photon polarization because the transversality condition behaves like a superselection rule. 
We can however define an effective reduced density matrix which corresponds to a restricted class 
of positive operator-valued measures. There are no pure photon qubits, and no exactly orthogonal 
qubit states. Reduced density matrices for the spin of massive particles are well-defined, but are 
not covariant under Lorentz transformations. The spin entropy is not a relativistic scalar and has 
no invariant meaning. The distinguishability of quantum signals and their entanglement depend on 
the relative motion of observers. 



I. INTRODUCTION 



The relationship between information and physics has 
been an intriguing problem for many years It took 
a new twist with the emergence of quantum informa- 
tion theory whose paradigm "information is physical" 
added a new point of view to old questions. Quan- 
tum information theory usually involves only a nonrela- 
tivistic quantum mechanics. However, both for the sake 
of logical completeness and in order to derive physical 
bounds on information transfer, its processing, and the 
errors involved, a full relativistic treatment is required. 
Additional motivation to relativistic extensions of quan- 
tum information theory comes from quantum cosmology. 
Quantum field theory in curved spacetime, and black hole 
physics in particular, present challenges that everybody 
who upholds the principle that "information is physi- 
cal" should respond to. Techniques of quantum infor- 
mation theorycan be applied to field theory and black 
hole physics In this presentation we describe some 
of the new features of quantum information theory when 
the effects of special relativity are taken into account. 

The concept of reduced density matrix is fundamental 
for quantum information. Its properties are significantly 
modified when we deal with relativistic effects. Let us use 
Latin indices for the description of a subsystem which is 
excluded from our description, and Greek indices for the 
subsystem that we can actually describe. The compo- 
nents of a state vector would thus be written Vm^i and 
those of a density matrix Pmfj,.ni'- The reduced density 
matrix of the system of interest is given by 

^^ly — ^ ^ Pm^.TLU • (1) 

m,n 

Even if /9 is a pure state (a matrix of rank one) , r is 
general a mixed state. Its entropy is defined as 

S = -tr(Tlogr). 



m 



(2) 



An important consequence of relativity is that there 
is a hierarchy of dynamical variables: primary variables 
have relativistic transformation laws that depend only 
on the Lorentz transformation matrix A that acts on 



the spacetime coordinates. For example, momentum 
components are primary variables. On the other hand, 
secondary variables such as spin and polarization have 
transformation laws that depend not only on A, but also 
on the momentum of the particle. As a consequence, the 
reduced density matrix for secondary variables, which 
may be well defined in any coordinate system, has no 
transformation law relating its components in different 
Lorentz frames. A simple example will be given below. 

Moreover, an unambiguous definition of the reduced 
density matrix by means of Eq. is possible only if 
the secondary degrees of freedom are unconstrained. For 
gauge field theories, that equation may be meaningless 
if it conflicts with constraints imposed on the physical 
states 0, 0- In the absence of a general prescription, 
a case-by-case treatment is required. A particular con- 
struction, valid with respect to a certain class of tests, 
is given in Sec. IIIII A general way of defining reduced 
density matrices for physical states in gauge theories is 
an open problem. 

The next two sections present detailed calculations 
and explore implications for distinguishability of quan- 
tum states and their entanglement. In the last section we 
show that a description of quantum channels by means of 
completely positive maps ji,6j is only an approximation. 



II. MASSIVE PARTICLES 



We first consider the relativistic properties of the spin 
entropy for a single, free particle of spin i and mass m > 
0. For massive particles a reduced density matrix is well- 
defined, but it has no invariant meaning 7]. The reason is 
that under a Lorentz boost, the spin undergoes a Wigner 
rotation |^ |^ 0| whose direction and magnitude depend 
on the momentum of the particle. Even if the initial 
state is a direct product of a function of momentum and 
a function of spin, the transformed state is not a direct 
product. Spin and momentum become entangled. 

The quantum state of a spin-i particle can be written, 
in the momentum representation, as a two-component 
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spmor, 



V'(P) 



ai(p) 
a2(p) 



(3) 



where the amphtudes satisfy J2r I \'^r{p)\'^dp = 1. 
The normaUzation of these amphtudes is a matter of con- 
venience, depending on whether we prefer to include a 
factor pq = (m^ + p^y^^ in it, or to have such factors in 
the transformation law © below. Here we shall use the 
second alternative, because it is closer to the nonrela- 
tivistic notation which appears in the usual definition of 
entropy. 

We emphasize that we consider normalizable states, 
in the momentum representation, not momentum eigen- 
states as usual in textbooks on particle physics. The lat- 
ter are chiefly concerned with the computation of (injout) 
matrix elements needed to obtain cross sections and other 
asymptotic properties. However, in general a particle has 
no definite momentum. For example, if an electron is 
elastically scattered by some target, the electron state 
after the scattering is a superposition that involves mo- 
menta in all directions. In that case, it still is formally 
possible to ask, in any Lorentz frame, what is the value of 
a spin component in a given direction (this is a legitimate 
Hermitian operator). 

Let us define a reduced density matrix, r = 
/ dp'i/;(p)V'^(p)j giving statistical predictions for the re- 
sults of measurements of spin components by an ideal 
apparatus which is not affected by the momentum of the 
particle. The spin entropy is 



S ^ -tr(rlogr) 



^A.log 



A, 



(4) 



where Xj are the eigenvalues of r. 

As usual, ignoring some degrees of freedom leaves the 
others in a mixed state. What is not obvious is that 
in the present case the amount of mixing depends on 
the Lorentz frame used by the observer. Indeed consider 
another observer (Bob) who moves with a constant ve- 
locity with respect to Alice who prepared state |(2Jl. In 
the Lorentz frame where Bob is at rest, the same spin-^ 
particle has a state 



(5) 



The transformation law is 



a'(p) = [(A-V)o/Po]'/' ^i?„[A,(A-V)]«.(A-V), 

S 

(6) 

where D^s is the Wigner rotation matrix for a Lorentz 
transformation A |3, 0, E| ■ 

As an example, take a particle prepared by Alice with 
spin in the z direction, so that 02 (p) = 0, and 



ai(p) = A^exp(-pV2A2), 



(7) 



where iV" is a normalization factor. Spin and momentum 
are not entangled, and the spin entropy is zero. When 



that particle is described in Bob's Lorentz frame, moving 
with velocity /3 in a direction at an angle 9 with Alice's z- 
axis, a detailed calculation shows that both a[ and a'2 are 
nonzero, so that the spin entropy is positive This 
henomenon is illustrated in Fig. 1. It can be shown 
[T]] | that a relevant parameter, apart from the angle 9, 
is in the leading order in momentum spread, 



m 



(8) 



where A is the momentum spread in Alice's frame. The 
entropy has no invariant meaning^ because the reduced 
density matrix r has no covariant transformation law, 
except in the limiting case of sharp momenta. Only the 
complete density matrix transforms covariantly. 




FIG. 1: Dependence of the spin entropy S, in Bob's frame, 
on the values of the angle 6 and a parameter T = [1 — (1 — 
/32)'/2]A/m/3. 

It is noteworthy that a similar situation arises for a 
classical system whose state is given in any Lorentz frame 
by a Liouville function . Recall that a Liouville func- 
tion expresses our probabilistic description of a classical 
system — what we can predict before we perform an 
actual observation — just as a quantum state is a math- 
ematical expression used for computing probabilities of 
events. 

Consider now a pair of orthogonal states that were 
prepared by Alice. How well can moving Bob distinguish 
them? We shall use the simplest criterion, namely the 
probability of error Pe, defined as follows: an observer 
receives a single copy of one of the two known states and 
performs any operation permitted by quantum theory in 
order to decide which state was supplied. The probability 
of a wrong answer for an optimal measurement is jl3j | 



PE{pi,P2) = \-\tr^{pi- P2?. (9) 

In Alice's frame Pe = 0. It can be shown that in 
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Bob's frame, oc F^, where the proportionaUty factor 
depends on the angle 6 defined above. 

An interesting problem is the relativistic nature of 
quantum entanglement when there are several particles. 
For two particles, an invariant definition of the entangle- 
ment of their spins would be to compute it in the Lorentz 
"rest frame" where (X]p) = 0- However, this simple def- 
inition is not adequate when there are more than two 
particles, because there appears a problem of cluster de- 
composition: each subset of particles may have a different 
rest frame. This is a difficult problem, still awaiting for 
a solution. We shall mention only a few partial results. 

Alsing and Milburn Q considered bipartite states 
with well-defined momenta. They showed that while 
Lorentz transformations change the appearance of the 
state in different inertial frames and the spin directions 
are Wigner rotated, the amount of entanglement remains 
intact. The reason is that Lorentz boosts do not cre- 
ate spin-momentum entanglement when acting on eigen- 
states of momentum, and the transformations on the pair 
are implemented on both particles as local unitary trans- 
formations which are known to preserve the entangle- 
ment. The same conclusion is also valid for photon pairs. 

However, realistic situations involve wave packets. For 
example, a general spin-i two-particle state may be writ- 
ten as 



CTl,(T2 



d/i(pi)d^(p2)g(CTlO-2, Pi, P2)|Pl, cri)«i|p2, ^2) 



where 



d^l{p) 



(10) 



(11) 



(27r)32pO' 

is a Lorentz-invariant measure. For particles with well 
defined momenta, g is sharply peaked at some values pio, 
P20- Again, a boost to any Lorentz frame S' will result 
in a unitary U{A) (g) U{A), acting on each particle sep- 
arately, thus preserving the entanglement. However, if 
the momenta are not sharp, so that the spin-momentum 
entanglement is frame dependent, then the spin-spin en- 
tanglement is frame-dependent as well. 

Gingrich and Adami 15] investigated the reduced den- 
sity matrix for IT12) and made explicit calculations for 
the case where g is a Gaussian, as in [3. They showed 
that if two particles are maximally entangled in a com- 
mon (approximate) rest frame (Alice's frame), then the 
degree of entanglement, as seen by a Lorentz-boosted 
Bob, decreases when the boost parameter /3 — > 1. Of 
course, the inverse transformation from Bob to Alice will 
increase the entanglement. Thus, we see that that spin- 
spin entanglement is not a Lorentz invariant quantity, 
exactly as spin entropy is not a Lorentz scalar. 



III. PHOTONS 

Relativistic effects that we describe in this section are 
essentially different from those for massive particles that 



were discussed above, because photons have only two 
linearly independent polarization states. The properties 
that we discuss are kinematical, not dynamical. At the 
statistical level, it is not even necessary to involve quan- 
tum electrodynamics. Most formulae can be derived by 
elementary classical methods Q . It is only when we con- 
sider individual photons, for cryptographic applications, 
that quantum theory becomes essential. The diffrac- 
tion effects mentioned above lead to superselection rules 
which make it impossible to define a reduced density ma- 
trix for polarization. As shown below, it is still possible 
to have "effective" density matrices; however, the latter 
depend not only on the preparation process, but also on 
the method of detection that is used by the observer. 

In applications to secure communication, the ideal sce- 
nario is that isolated photons (one particle Fock states) 
are emitted. In a more realistic setup, the transmission 
is by means of weak coherent pulses containing on the 
average less than one photon each. A basis of the one- 
photon space is spanned by states of definite momentum 
and helicity, 

(12) 



|k,4)^|k) 



-k /' 



where the momentum basis is normalized by (q|k) = 
(27r)3(2/e°)(5(3)(q - k), and helicity states |e^) are ex- 
plicitly defined by Eq. H16|) below. 

As we know, polarization is a secondary variable: 
states that correspond to different momenta belong to 
distinct Hilbert spaces and cannot be superposed (an ex- 
pression such as je^ ) + l^q ) meaningless if k 7^ q). The 
complete basis ()12|l does not violate this superselection 
rule, owing to the othogonality of the momentum basis. 
Therefore, a generic one-photon state is given by a wave 
packet 



1^1-)= I d^i{k)f{k)\k,a{k)). 



(13) 



The Lorentz-invariant measure is dfi{k) = d'^'k/ (2Tr)^2k'^ , 
and normalized states satisfy / d/i(A:)|/(k)p = 1. The 
generic polarization state |Q!(k)) corresponds to the geo- 
metrical 3-vector 



Q;(k) = a^{k)e^ + Q;_(k)e]^ . 



(14) 



where |q;+P + — 1, and the explicit form of Cj, is 

given below. 

Lorentz transformations of quantum states are most 
easily computed by referring to some standard momen- 
tum, which for photons is p''' — (1,0,0, 1). Accordingly, 
standard right and left circular polarization vectors are 
— (1, ±i, 0)/\/2. For linear polarization, we take 
Eq. 114|) with a+ = («-)*, so that the 3- vectors a.{k) 
are real. In general, complex cy.{k) correspond to elliptic 
polarization. 

Under a Lorentz transformation A, these states become 
|kA, a(kA)), where kA is the spatial part of a four- vector 
^A = Afc, and the new polarization vector can be ob- 
tained by an appropriate rotation |14L Il6| 



a(kA) = i?(kA)i?(k)-'a(k), 



(15) 
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where k is the unit vector in the direction of k. Finally, 
for each k a polarization basis is labeled by the helicity 
vectors, 



4 = RO^H 



(16) 



Let us try to define a reduced density matrix in the 
usual way, 

p = yd/i(k)|/(fc)p|k,a(k))(k,a(k)|? (17) 

The superselection rule that was mentioned above does 
not forbid this definition, because only terms with the 
same momentum k are summed. However, since polar- 
ization is a secondary variable, this object cannot have 
definite transformation properties under boosts. This de- 
ficiency is familiar to us from the analysis of reduced den- 
sity matrices of massive particles. However, for mass- 
less particles, the situation is worse: POVMs that are 
given by 2 X 2 matrices represent measurement devices 
and should transform under a representation of the rota- 
tion group 0(3). On the other hand, even for the com- 
plete photon state, ordinary rotations of the reference 
frame correspond to elements of E(2), so that probabili- 
ties would not be invariant under rotations. 

Therefore, let us find a more physical definition of a re- 
duced density matrix for polarization J^. The labelling 
of polarization states by Euclidean vectors ej^ suggests 
the use of a 3 x 3 matrix with entries labelled x, y and 
z. Classically, they correspond to different directions of 
the electric field. For example, a reduced density matrix 
Px would give the expectation values of operators rep- 
resenting the polarization in the x direction, seemingly 
irrespective of the particle's momentum. 

To have a momentum-independent polarization is to 
tacitly admit longitudinal photons. Unphysical concepts 
are often used in intermediate steps in theoretical physics. 
Momentum-independent polarization states thus consist 
of physical (transversal) and unphysical (longitudinal) 
parts, the latter corresponding to a polarization vector 
= k. For example, a generalized polarization state 
along the a;-axis is 



|x)=x+(k)|e+; 



-(k)|* 



a;£(k)|< 



k/i 



(18) 



where a;±(k) = • x, and ^^(k) = x • k = sin cos 0. It 
follows that |a;_|_ 



^ ' ja^-P -I- \xi\'^ = 1, and we thus define 



ex{k) = 



2;+(k)e^ +a;_(k)ei^ 



(19) 



as the polarization vector associated with the x direction. 
It follows from (fTH|l that (x|x) = 1 and (x|y) = x-y = 0, 
and likewise for other directions, so that 

|x)(i| + |y)(y| + |z)(z|==lp, (20) 

where Ip is the unit operator in momentum space. 



We can now define an "effective" reduced density ma- 
trix adapted to a specific method of measuring polar- 
ization, as follows Q. To the direction x corresponds a 
projection operator 

= |X) (5c| §5 Ip = |i) (i| ® J dfi{k)\k) (k| , (21) 

The action of on j^*) follows from Eq. ((TH|l and 
(^kl^k) ~ 0- Only the transversal part of |x) appears 



in the expectation value: 



(vI^|P,|vI')= J dp{k)\f{k)f\x+{k)a*,{k)+x-{k)a*l{k)\' 



(22) 



It is convenient to write the transversal part of |x) as 



|b.(k)) ^ (I 



-k/\^k I ^ i"=k/\"=k 



ekX^kDIx) 



x+(k)|e+)+x_(k)|ek). 



(23) 



Likewise define |bj,(k)) and |b2(k)). These three state 
vectors are neither of unit length nor mutually orthogo- 
nal. For k = (sin 9 cos (p, sin 9 sin ip, cos 9) we have 

|b,(k)) =c(0,0)|k,e,(k)), (24) 

where e^ik) is given by Eq. H19|l . and c{9,(j)) = 
^xl+xl. 

Finally, a POVM element which is the physical 
part of Px , namely is equivalent to Px for physical states 
(without longitudinal photons) is 



Ex = J d^l(k)\k,hx{k)){k,hx{k)\, 



(25) 



and likewise for other directions. The operators Ex, Ey 
and E^ indeed form a POVM in the space of physical 
states, owing to Eq. H2U|) . The above derivation was, ad- 
mittedly, a rather circuitous route for obtaining a POVM 
for polarization. This is due to the fact that the latter 
is a secondary variable, subject to superselection rules. 
Unfortunately, this is the generic situation. 

To complete the construction of the density matrix, 
we introduce additional directions. Similarly to a stan- 
dard practice of channel matrices reconstruction [T^ , we 
consider Ex+z.x+z, Ex~iz,x-iz and similar combinations. 
For example. 



E. 



X'{-Z,X-\-Z 



K|x)-f |z))((x| + (z|)®lp 



(26) 



Let us denote |x)(z| (8) Ip as Exz, even though this is 
not a positive operator. We then get a simple expres- 
sion for the reduced density matrix corresponding to the 
polarization state |Q:(k)): 

Prnn = (*|£;,„„|*) = 

d/x ( A) I / (k) 1 2 (a (k) I (k)) (b„ (k) I a (k)) , (27) 
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where m,n,= x,y,z. It is interesting to note that this 
derivation gives a direct physical meaning to the naive 
definition of a reduced density matrix, 

Pml™ = / dA*(k)|/(fc)pa™(k)a:(k) = p^n (28) 

Our basis states |k, et) are direct products of momen- 
tum and polarization. Owing to the transversality re- 
quirement ek • k = 0, they remain direct products un- 
der Lorentz transformations. All the other states have 
their polarization and momentum degrees of freedom en- 
tangled. As a result, if one is restricted to polarization 
measurements as described by the above POVM, there do 
not exist two orthogonal polarization states. An immedi- 
ate corollary is that photon polarization states cannot be 
cloned perfectly, because the no-cloning theorem |0 for- 
bids an exact copying of unknown non-orthogonal states. 
In general, any measurement procedure with finite mo- 
mentum sensitivity will lead to the errors in identifica- 
tion. First we present some general considerations and 
then illustrate them with a simple example. 

Let us take the z-axis to coincide with the average di- 
rection of propagation so that the mean photon momen- 
tum is fc^z. Typically, the spread in momentum is small, 
but not necessarily equal in all directions. Usually the 
intensity profile of laser beams has cylindrical symmetry, 
and we may assume that A^^ ~ ~ A^ where the index 
r means radial. We may also assume that A^^A^. We 
then have 

/(k) CX /i[(fc, - fc^)/A,] f2{kr/Ar). (29) 

We approximate 

9 w tanff = kr/kz ~ kr/kA- (30) 

In pictorial language, polarization planes for different 
momenta are tilted by angles up to ~ A^/fc^, so that 
we expect an error probability of the order A^/fc^. In 
the density matrix pmn sll the elements of the form pmz 
should vanish when A^ — > 0. Moreover, if A^ — > 0, 
the non-vanishing xy block goes to the usual (monochro- 
matic) polarization density matrix. 

As an example we consider two states which, if the mo- 
mentum spread could be ignored, would be |A:^z,e^^.). 
To simplify the calculations we assume a Gaussian dis- 
tribution: 

/(k) = iVe-(^-^-^--^)'/2A^e-^-'/2^', (31) 

where iV is a normalization factor and A^ <C A^. More- 
over, we take the polarization components to be = 

R(k)e^. That means we have to analyze the states 

|*±)=y'dMk)/(k)|e± k), (32) 
where /(k) is given above. 



It can be shown |5| that in the leading order 

PE{p+,P-)^Al/Ak\. (33) 

Now we turn to the distinguishability problem from the 
point of view of a moving observer. Bob. The probability 
of an error by Bob is still given by Eq. 0. The distin- 
guishability of polarization density matrices depends on 
the observer's motion. We again assume that Bob moves 
along the z-axis with a velocity v. Detailed calculations 
m show that 

P'e = \^^Pe, (34) 

which may be either larger or smaller than Pe- As ex- 
pected, we obtain for one-photon states the same Doppler 
effect as in the classical calculations 0. 

Since maximally entangled states are pure, the fact 
that all polarization density matrices are mixed implies 
that maximal EPR-type correlations shall never be ob- 
served, and that maximal attainable value will depend 
on the momentum spread of the states. 

IV. COMMUNICATION CHANNELS 

Although reduced polarization density matrices have 
no general transformation rule, the above results show 
that such rules can be established for particular classes of 
experimental procedures. We can then ask how these ef- 
fective transformation rules, t' = T{t), fit into the frame- 
work of general state transformations. Equation H35|) 
gives an example of such a transformation. 

p' = p{l- —) + {a^p(7^ + aypay) — . (35) 

It relates the reduced density matrix p, obtained from 
Alice's Eq. to Bob's density matrix p' . This partic- 
ular transformation is completely positive; however, it is 
not so in general. 

It can be proved that distinguishability, as expressed 
by natural measures like Pe , cannot be improved by any 
completely positive transformation [l3l | . It is also known 
that the complete positivity requirement may fail if there 
is a prior entanglement with another system |l9j | . Since in 
the two previous sections we have seen that distinguisha- 
bility can be improved, we conclude that these transfor- 
mations are not completely positive. The reason is that 
the Lorentz transformation acts not only on the "inter- 
esting" discrete variables, but also on the "hidden" mo- 
mentum variables that we elected to ignore and to trace 
out, and its action on the interesting degrees of freedom 
depends on the hidden ones. 

This technicality has one important consequence. The 
notion of a channel is fundamental both in classical and 
quantum communication theory. Quantum channels are 
described as completely positive maps that act on qubit 
states 20]. Qubits themselves are realized as particles' 
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discrete degrees of freedom. If relativistic motion is im- quantum channel, but the very representation of a chan- 
portant, then not only does the vacuum behave as a noisy ncl by a completely positive map fails jj]. 



[1] E. T. Jaynes, Phys. Rev. 108, 171 (1958); J. A. Wheeler 
and W. H. Zurek, editors. Quantum Theory and Mea- 
surement (Princeton University Press, Princeton, 1983). 

[2] R. Landauer, Physics Today 44 (5), 23 (1991). 

[3] A. Peres and D. R. Terno, Rev. Mod. Phys. 75 (2003), 
in press; e-print quant-ph/0212023 

[4] D. Beckman, D. Gottesman, A. Kitaev, and J. Preskill, 
Phys. Rev. D 65, 065022 (2002). 

[5] A. Peres and D. R. Terno, J. Mod. Opt. 50 (2003), in 
press, e-print quant-ph/0208128 

[6] E. B. Dsy'ies, Quantum Dynamics of Open Systems (Aca- 
demic Press, New York, 1976); K. Kraus, States, Effects, 
and Operations: Fundamental Notions of Quantum The- 
ory (Springer, Berlin, 1983). 

[7] A. Peres, P. F. Scudo, and D. R. Terno, Phys. Rev. Lett. 
88, 230402 (2002). 

[8] S. Weinberg, The Quantum Theory of Fields (Cambridge 
University Press, Cambridge, 1995) vol. 1. 

[9] F. R. Halpern Special Relativity and Quantum Mechanics 
(Prentice-Hall, Englewood Cliffs, 1968). 
[10] N. N. Bogolubov, A. A. Logunov, A. I. Oksak, and I. T. 
Todorov, General Principles of Quantum Field Theory 
(Kluwer, Dordrecht, 1990). 



[11] D. R. Terno, Quantum Information and Relativity The- 
ory, Ph.D. thesis (Technion, Haifa, 2003). 

[12] R. Balescu and T. Kotera, Physica (Utrecht) 33, 558 
(1967). 

[13] C. A. Fuchs and J. van de Graaf, IEEE Trans. Info. The- 
ory IT-45, 1216 (1999). 

[14] P. M. Alsing and G. J. Milburn, Quant. Info. Comp. 2, 
487 (2002). 

[15] R. M. Gingrich and C. Adami, Phys. Rev. Lett. 89, 
270402 (2002). 

[16] D. Han, Y. S. Kim, and D. Son, Phys. Rev. D 31, 328 
(1985). 

[17] I. L. Chuang and M. A. Nielsen, J. Mod. Opt. 44, 2455 
(1997). 

[18] W. K. Wootters and W. H. Zurek, Nature 299, 802 
(1982); D. Dieks, Phys. Lett. A 92, 271 (1982). 

[19] P. Pechukas, Phys. Rev. Lett. 73, 1060 (1994); P. 
Stelmachovic and V. Buzek, Phys. Rev. A 64, 062106 
(2001). 

[20] A. S. Holevo, Russ. Math. Surveys 53, 1295 (1999); C. 
King and M. B. Ruskai, IEEE Trans. Info. Theory IT-47, 
192 (2001); M. Keyl, Phys. Rep. 369, 431 (2002). 



